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ABSTRACT 



We construct a structural model of the Andromeda Galaxy, simultaneously corre- 
sponding to observed photometrical and kinematical data and chemical abundances. 
In this paper we present the observed surface brightness, col our and metallicity distri- 
butio ns, and compare them to the model galaxy. In Paper II (|Tempel. Tamm &: Teniei 
120071) we present similar data for the kinematics, and derive the mass distribution of 
the galaxy. 

To determine structural parameters of stellar components, we have collected the 
observed U, B, V, R, I and L luminosity distributions and have calculated the re- 
lated colour indices. By using far-infrared imaging data of M 31 and a thin dust disc 
assumption, we derive dust-free surface brightness and colour distributions. 

The model galaxy is constructed as a superposition of four axially symmetric stel- 
lar components: a bulge, a disc, an inner halo and an extended diffuse halo. A set of 
colour indices and metallicity is ascribed to each component, consistent with chemical 
evolution models of simple stellar populations, metallicity distribution measurements 
and star formation history estimates. The models of chemical evolution allow to esti- 
mate mass-to-light ratios of the stellar components. 

We find the total absorption corrected luminosity of M 31 to be Lb — (3.3 ± 0.7) • 
10^°Lq, corresponding to an absolute luminosity Mb = —20.8 ± 0.2 mag. Of the total 
luminosity, 41 % (0.57 mag) is obscured from us by the dust inside M31. The intrinsic 
colours of the stellar populations are in agreement with an extended formation time 
of bulge and disc stars 6-10 gigayears ago, with an additional population of young 
stars making up about 10 % of the disc luminosity and less than 1 % of its mass. Halo 
colours are not well reproduced by chemical evolution models; especially in the outer 
halo, the actual colours are considerably redder than the modelled ones. The chemical 
evolution model yields the following mass-to-hght ratios: M/Lb =3.7-6.9 Mq/Lq for 
the bulge, M/Lb =2.9-5.2 M0/L0 for the disc, M/Lb =2.7-5.7 Mq/'Lq for the 
inner halo and M/Lb ~ 2 Mq/Lq for the outer diffuse halo. The total mass of the 
visible matter is M^is =(10-19)T0"'^°Mq, giving the total intrinsic mass-to-light ratio 
of the visible matter M/Lb =3.1-5.8Mq/Lq. 

The use of the model parameters for a dynamical analysis and for determining 
dark matter distribution is presented in Paper II. 

Key v^rords: galaxies: individual: Andromeda, M 31 - galaxies: photometry - galaxies: 
stellar content - galaxies: structure - dark matter 



1 INTRODUCTION 

Models of hierarchical merging in a cold dark matter 
(CDM) universe describe rather well the observed large- 
scale structure of clusters, superclusters and their net- 
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work. They are also successful in reproducing general 
aroperties of galaxies with different morphological types 



properties or galaxies with dinerent morphological types 
jAbadi et all l2003l: iBell et al.l l2003l : iNagamine et al.l l2005l : 
iGovernato et al.ll2007f K 



On the other hand, the observed number of dwarf 
galaxies seems to be too small in comparison with CDM 
models l|Klvpin et al.ll 19991 : [Moore et al.|[l999l : iBenson et all 
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|2002| ). In the case of the Ursa Minor dwarf spheroidal galaxy, 
the observed velocity dispersions are in contradiction to th e 
cuspy density profile of CDM haloes (jKlevna et al.ll2003h . 
But by far the most commonly referred contradiction is 
related to central densities of dark matter halos. Cosmolog- 
ical simulations generate DM hal oes with central density 
cusps of p ^ o r steep er (e.g. iNavarroT Frenk fc Wh"it3 
1 19971 : iMoore et all Il999l ). Observations of dwarf and 
low surface brightness disc galaxies have usually shown 



than cuspy profiles (|Blais-Ouellette. Amram & Carignan 


2OOII: iBorriello & Salucci 200ll: de Blok & Bosma 


2OO2I: ISalucci. Walter & Borrielld 


2003; 


Weldrake. de Blok & Walter 20031: Simon et al.l 


200a 


Zackrisson et al.l 


20061: iKassin. de Jone & Weined 


2009: 


Gentile et al. 2004 


.l2007l:IValenzuela et al.ll2007f). 



For a proper analysis of actual galaxies, it is necessary 
to know the distribution of both visible and dark matter. 
Unfortunately, the structure and mass distribution of stel- 
lar populations is not known precisely enough. One of the 
most uncertain aspects here is the mass-to-light ratio (AI/L) 
of visible matter. In principle, chemical evolution models al- 
low to calculate the evolution of both luminosity and mass 
of a stellar population, but they involve several insufficiently 
constrained parameters (the Initial Mass Function, starting 
time, duration and intensity of star formation in difi'erent 
subsystems, initial metallicities etc). It is possible to de- 
crease the degeneracy in M/L determinations by using 3-5 
different colour indices, but due to measurement uncertain- 
ties, colour indices are often controversial. Internal dust at- 
tenuation affects the surface brightness distribution profiles 
and colour indices, causing additional uncertainty in the dis- 
tribution of visible matter. 

A comparison of the results of galactic structure mod- 
elling with additional and, moreover, independent observa- 
tions (metallicities of stellar populations, stellar rotation 
curve, velocity dispersions along several slit-positions) al- 
lows us to constrain the distributions of visible matter and 
thereafter DM. 

In this work, we analyse the density distributions of vis- 
ible and DM components in a nearby luminous disc galaxy, 
the Andromeda galaxy M31. The galaxy M31 was selected 
because (1) photometrical and kinematical (rotation, dis- 
persions) data are known with sufficiently high resolution 
in order to study the bulge region; (2)velocity dispersions 
have been measured also outside the galactic apparent ma- 
jor ax;is; (3) direct measurements of stellar metallicities al- 
low to constrain M/L of visible matter; (4) independent 
estimates of mass distribution on large scales are avail- 
able (globular clusters, satellites, stream, kinematics of the 
Milky Way -|- M 31 syst em) ; (5) in the case of the Milky Way, 
iBinnev fc Evand HooJ) derived that cuspy CDM haloes are 
inconsist ent with the obs e rvatio nal data (but see a different 
study bv lBattagha et al.l (120051 )). 

Composite models of galaxies (including M31), taking 
into account surface photometry and the history of chem- 
ical evolutio n and kinematical data, were first constructed 
by Einast^ (119741 ). It was demonstrated that these models 
allow to distinguish stellar populations in galaxies and to 
calculate their main parameters. In addition to this kind of 
general mass and luminosity distribution models, t he struc- 
ture of M 31 stellar populations have been studied bv lHoddeil 



(1 19951). who adopted iBahcall fc Soneiral (| 19841 ) model, and 
bv lBellazzini et all (|2003l ). A rather different kind of model 
on the basi s of phase-space d istribution function was con- 
structed bv lWidrow fc Dubinsk i (2005). 

In this paper, we construct a photometric model of M 31 
stellar populations on the basis of s urface brightness profiles 
in U, B, V, R, I, L colours (see also lTenies. Hand fc Einastol 
I1994'). The derived photometrical model gives us parameters 
of galactic components, colour indices among them. Stellar 
metallicities are available from independent spectral obser- 
vations and colour-magnitude diagrams. We have used all 
these data (component radii, luminosities, colour indices, 
metallicities) as input parameters for chemical evolution 
models to calculate the ages and M/Ls of the components. 
The main output of the present paper are density distribu- 
tion parameters and M/Ls of visible g alactic c omponents. 

In the companion paper (Paper II, Temp el et al. 2007i), 
we apply the results of the present paper and construct a 
mass distribution model of M 31. Calculating the rotation 
velocities and velocity dispersions o f visible matter with the 
help of the dynamical model (see iTempel fc Tenie"3 12006| ) 
we can estimate the amount of DM which must be added 
to obtain an agreement with the measured rotation and dis- 
persion data. 

We have applied the following general param- 
eters in our calculations: the inclination ang le of 

M31 has been taken 77.5 ° (jWalterbos fc KennicuttI Il988l : 
Ide Vaucouleurs et al.lll99l | ). the major axis position angle 
is 38.1° (|Walterbos fc Kennicutt|[l987] : iFerguson et al.ll2002l) 
and the distance to M31 is 785 kpc ( McConnachie et al.l 
l2005h . corresponding to the scale 1 arcmin = 228 pc. Ab- 
sorption in the Milky Way has b een taken according to 
ISchlegel. Finkbeiner fc David (|1998D . 



2 OBSERVATIONAL DATA 

In this Section, we describe observational data used to con- 
struct a photometrical model and to determine metallicities 
of galactic components. 



2.1 Photometrical data 

By now surface photometry of M31 is available in U, B, 
V , R, I , J, H , K , L colours. Referen ces to earlier observa- 
tions ar e given in iTenies et al.1 (Il994l). In outer parts of the 
galaxy, IPritchet fc van den Bergh (l99j ) and llrwin et al.l 
(20051) derived surface brightness profiles along the minor 
axis. The K-profile measu r ed on the basis of star counts by 
IPritchet fc van den Berghl (|l994l ) extends up to the surface 
brightness 29.3 mag arcsec"^. They estimated a l so th e axial 
ratio of isophotes in outer regions, llrwin et al.1 (|2005l ) mea- 
sured surface brightness in V and / colours directly and 
with the help of star counts and compiled profiles up to 
31.3 mag arcsec"^ in V and up to 29.2 mag arcsec"^ in /. 

Recent observat ions with the Spitzer Space Telescope 
(|Barmbv et al.ll2006l ) have allowed to extend the photome- 
try to the i-passband. The J, H , Jf-profiles have a rather 
limited spatial extent and resolution, and we decided not to 
use them in here. The wavelengt h of the Galaxy Ev olution 
Explorer near- UV observations (|Thilker et al.ll2005l ) is un- 
fortunately too much off the standard U for interpretation 
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Upper lines: extinction corrected 
Lower lines: original observations 
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Figure 1. M31 luminosity distribution along the major semiaxis 
in , B, V, R, I, L colours. For each colour, the lower distribu- 
tion corresponds to the measured data points (Section l2.lt . and 
the upper distribution presents our surface brightness distribu- 
tion corrected from the intrinsic absorption (Section 13.11 1. Note 
that except for U, all other distributions have been shifted along 
the j/-axis as indicated in the legend. 



with chemical evolution models and we have not used these 
data here. 

The composite surface brightness profiles in U, B, V, 
R, I, L colours along the major and/or the minor axes were 
derived by averaging the results of different authors. Dif- 
ferent R and / colour system profiles were tr ansferred into 
the C ousins system, using the calibration by iFrei fc GunnI 
l|l994 ). The resultant observational profiles along the ma- 
jor axis together with the estimated errors are presented in 
Fig. [1] All the surface brightness profiles obtained in this 
way belong to the initial data set of our model construction. 



2.2 Metallicity measurements 

To constrain the chemical evolution models, we have col- 
lected from the literature most recent metallicity estimates 
at various locations over the galaxy, extending from the cen- 
tral bulge region to the outermost halo fields; the sources are 
listed in the legend of Fig. [21 Within these works, the metal- 
licities of the dense bulge and disc regions have been esti- 
mated by fitting isochrones to colour-magnitude diagrams of 
the red giant branch (RGB) stars found in the fields. In the 
outer low-density fields, isolating of M 31 stars from the fore- 
ground Milky Way dwarfs becomes a crucial issue, and spec- 
troscopic measurements are necessary for higher reliability. 
Furthermore, according to recent studies, the outer regions 
of M31 contain a number of distinct stellar populations in 
the form of stellar streams, which probably are remnants of 
for mer satellites of M 31 with their own distinct metallic- 
ity (jlbata et al.ll2007l V Thus metallicity estimates based on 
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Figure 2. Metallicities of M 31 along the major axis and the mi- 
nor axis. Open symbols represent photometric estimates, filled 
symbols represent spectroscopic measurements. The solid line 
gives the best-fit metallicity distribution of our model; the shaded 
region shows the whole possible metallicity range, used to deter- 
mine the uncertainties of the modelled colours and mass-to-light 
ratios (Section 14.11 1. 



randomly selected outer fields may introduce a significant 
scatter and even a bias. 

In difi'erent papers, metallicity has been presented in 
both iron-to- hydrogen ratio [Fe/H] as well as metal-to- 
hydrogen ratio [M/H] (also called the "global metallic- 
ity" scale). Photometric properties of a star depend on 
the abundances of all heavy elements rather than iron 
alone, thus we have converted all the measurements into 
the [M/H] scale as a more informative one. Conversion 
from one scale to another is not straightforward; for ex- 
ample, in the Milky Way disc, solar-metallicity stars have 
been shown to have metal-to-iron ratio close to that 
of the Sun, while an enhancement of the abundance of 
the a-elements with respect to iron occurs with a de- 
creasing metallicity an d thus with an increasing radius 
l|Edvardsson et all 1 19931 ): a similar enhancement has been 
found for the Mil ky Way halo stars (e.g. IZhao fc MagainI 
199(1). Following [Salaris, Chieffi fc Straniero l (119931 ) aiid 
Ferraro. Valenti fc Origlial ( 2006il ). in the disc region we ap- 
ply a conversion [M/H] = [Fe/H] + log(0.638/a + 0.362), 
where the a-element enhancement factor fa = 10^" '^''^'^''^' 
within the range > [Fe/H] > —1. In the halo fields, we 
apply [M/H] = [Fe/H] -I- 0.2, assuming the a-element en- 
hancement by a factor of 2. 

In Fig. [2l open symbols represent photometric metallic- 
ity estimates and filled symbols stand for spectroscopic mea- 
surements. The data points in the upper and lower panel of 
Fig. [Hare given as projected along the major and minor axis 
of the galaxy, respectively, accounting for the mean elliptic- 
ity of the galaxy at the location of each metallicity mea- 
surement. Due to differences in the geometry of the galactic 
components, such a projection cannot be correct for fields 
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located too much ofT the axes onto which the projection is 
made, thus we have not used data from disc fields deviat- 
ing from the target axis by angles more than about 20 de- 
grees. The error bars of metallicities present both the direct 
measurement uncertainties and variations arising from the 
uncertainty of t h e dist ance to the galajcy , redde ning etc; see 
iBellazzini et alj (120031 ) and iBrown et all (|2006h for an esti- 
mate of these effects. To consider all the uncertainties of our 
model resulting from this rather vaguely known metallicity 
distribution, we have used in our calculations the whole al- 
lowed range of metallicities, presented with the grey area in 
Fig.E 



3 SURFACE BRIGHTNESS DISTRIBUTION 
MODEL 

3.1 Recovering dust- free luminosity distributions 

Before fitting any models to the observations it is essen- 
tial to consider possible effects of dust extinction on the 
observed data. Ignoring of dust effects may lead to substan- 
tially wrong estimates of the structure, luminosity, colour 
and the mass of the stellar components. Accounting for the 
effects of intrinsic dust extinction, however, is a complicated 
task. Fortunately, the Andromeda is an intensely studied 
galaxy and its dust distribution has been revealed by sur- 
veys with IRAS, ISO and Spitzer satellites at a variety of 
far-infrared wavelengths. 

Conversion of the infrared flux along each line-of-sight 
into absorption is far from trivial, even if we assume the 
flux to be proportional to the dust column density, and we 
know the absorption properties of the dust. Also, the actual 
absorption along each line-of-sight strongly depends on the 
geometry of the stellar population s and the dust l ayer. The 
Spitzer telescope infrared view in ICordon et all ||2006D re- 
veals dust concentrated in several spiral segments in the disc 
region of the galaxy, peaking at about a dozen kiloparsecs 
from the centre. 

Relying on reddening measured f rom colour histogram s 
of globular clusters in M 31 halo by iBarmbv" et al.1 l|2000t ). 
we have assumed that the extinction law of M 31 is simi- 
lar to that of the Galactic dust. To calculate the total ab- 
sorption, we use optical depth estimates from the literature. 
IXu fc Heloul (|l996i ) constructed a sophisticated model to es- 
timate face-on optical depth of the M 31 dust on the basis of 
UV, optical and infrared luminosity, the latter relying on the 
IRAS measurements at 60 ^im and 100 \im. They found the 
azimuthally averaged optical depth tv to be in t he range of 
0.7-1. 6 at 2-14 kpc from the centre. Later on, iHaas et al.l 
added the ISO satellite measurements at 175 \im. 
By transferring the infrared emission into dust mass and 
the latter into optical depth, they r eache d an optical depth 
value close to that of lXu fc Heloul (|l996l ). However, to cor- 
rect optical luminosity profiles, the uncertainties are large, 
further increased by our limited knowledge about the 3-D 
distribution and dumpiness of the dust clouds. While in 
the disc region it is reasonable to consider dust being con- 
centrated within a rather thin dust disc as a result of disc 
kinematics, in the central bulge regions the dust clouds are 
likely to move along more varied orbits. A spherical rather 
than discy distribution of dust is suggested also by the far- 



infrared images. Unfortunately, accounting for the absorp- 
tion by roughly spheroidally distributed dust is too uncer- 
tain, thus we have avoided using the data within about 2 
kiloparsecs from the centre along the minor axis. 

The recovery of the major axis profiles turns out to be 
a much simpler task. We assumed the dust disc to be thin 
compared to the stellar disc, with its optical depth propor- 
tional to the infrared flux, and calculated the absorption 
along each line-of-sight. Interestingly, the resultant absorp- 
tion is rather insensitive to the assumptions about the dust 
disc thickness and face-on optical depth. The high inclina- 
tion angle of the galaxy increases the effective optical depth 
of the dust disc so much that nearly all the optical signal 
coming from the other side of the dust layer become ab- 
sorbed, and naturally all the light from the stars between us 
and the dust layer reaches us unaffected. Thus in the case 
of a thin dust layer with respect to the stellar disc, along 
the major axis absorption is effectively grey, being about 
0.75 mag almost uniformly for U, B, V, R, I filters and for 
any optical depth within the range tv = 0.5 — oo. The orig- 
inal and restored lum inosity p rofiles are presented in Fig. [T] 

For comparison, [Driver e t al. (2000) have derived in- 
tegral absorption value as a function of inclination angle 
statistically on the basis of over 10,000 galaxies. For the in- 
clination angle of M31, the absorption is 0.75-0.85 mag for 
the bulge and 0.6-0.65 mag for the disc in 5-filter, which 
matches well with our estimate of 0.75 mag. 

3.2 Choice of components 

A satisfactory fit to the measured luminosity distribution of 
M 31 can be achieved with just one or two Sersic profiles 
if only one-dimensional p hotometry, based o n a single filter 
observations, is used (e.g. Ilrwin et al.ll2005l ). In the case of 
M31, much more information is available: multi-colour pho- 
tometry along both major and minor axes, metallicity maps, 
kinematical data - a comprehensive model should attempt 
to fit all these observations. A better fit can be achieved 
by increasing the number of components, for example it has 
been revealed from the infrared maps of M 31 that the dusty 
star-forming region actually consists of several rings wit h 
offset centres and inclination angles (jCordon et al.l 120061 ). 
However, the use of a large number of components would 
blur our vision about the structure and the origin of the 
galaxy, thus here we try to use a minimal number of compo- 
nents giving simultaneously a satisfactory fit to all collected 
data. 

While the classical galactic components of the bulge and 
the disc are clearly visually detectable, and their existence 
and structure can undoubtedly be confirmed by studying 
the morphology of the galaxy and the kinematics of the em- 
bedded stars, classification of the outer stars into distinct 
stellar components is much more complicated. The outer re- 
gions are of extremely low density and the foreground stars 
of the Milky Way dominate in the observed fields. More- 
over, recent studies indicate the presence of a considerable 
non-relaxed substructure around M 31 in the form of tidal 
streams and disrupted companions (|lbata et al.ll2007l ). the 
most prominent of which is the Giant Steam. 

Several observational studies of distant fields along the 
minor axis have found evidence for the existence of stel- 
lar populations far beyond the bulge and the disc regions 
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l|lrwin et al.ll2005l : llbata et al.ll2007l '). The outermost regions 
seem to be dominated by a faint, diffuse, very extended 
metal p oor fialo, pe rfiaps embedding most of tfie wliole Local 
Group l|lbata et al. 2007). We refer to tliis underlying pop- 
ulation as the "outer diffuse halo". Between the radii from 
a few to about 25 kiloparsecs, recent studies suggest the 
presence of a dynamically warm but slowly rotating stellar 
population of stars with intermediate metallicity, which we 
designate as the "inner halo" . Note that t he same region has 
been referred to as th e "extended disc" (llbata et al]|2005l . 
I2OO7I'). the "spheroi d" (iBrown et al.ll2007h and the "bulge" 
i Kalirai et al.|[200^ ). depending on the viewpoint. Earlier, 



when no information about the underlying metal-poor pop- 
ulation was available, this population was simply referred to 
as the "halo". 

We have not attempted to include the Giant Stream in 
our model. The Giant Stream stars lie slightly off the minor 
axis of the galaxy and thus do not contribute significantly 
to our luminosity and kinematics data, while its expected 
total mass and luminosity can safely be fitted within the un- 
certainties of our results. Besides, our model can deal with 
axially symmetric stellar populations only. For similar rea- 
sons, we also neglect the suggested bar population and the 
elongated (double) nucleus of the galaxy. 

For simplicity, in our model we ascribe a set of a chem- 
ical composition and a formation time and subsequently a 
metallicity, M/Ls and colour indices homogeneously over 
each component. In reality, the parameters listed above are 
all, to some extent, radius-dependent. However, the usage of 
a superposition of several components enables to compen- 
sate for such a simplification and provides a fairly good fit 
to the actual observational data, as will be shown below and 
in Paper 11. 



3.3 Model density distribution of the components 

To construct a sufficiently flexible and consistent model, al- 
lowing to describe both the surface luminosity distribution 
of components and their dynamics, we start from a spatial 
density distribution law for individual components, which 
allows an easier fitting simultaneously for light distribution 
and kinematics. 

In such a model, the visible part of the galaxy is given 
as a superposition of its individual stellar components. The 
spatial density distribution of each visible component is 
approximated by an inhomogeneous ellipsoid of rotational 
symmetry with the constant axial ratio q and the density 
distribution law 



l{a) = 1(0) exp 



(—) 

\kaoJ 



l/N 



(1) 



where l{0) — hL/{4nqaQ) is the central density and L is the 
component luminosity; a = ^ B? -\- jq^, where R and z 
are two cylindrical coordinates; ao is the harmonic mean 
radius which characterizes rather well the real extent of 
the component, independently of the parameter A'^. The co- 
efficients h and k are normalizing parameters, depending 
on A'', which allows the density behaviour to vary with a. 
The definition of the normalizing parameters h and fc and 
their c alculation is described in Appendix B of lTenies et al.l 
l|l994 ) . The luminosity density distribution ([ij proposed in- 



dependently bv lEinastol (|l969l ) is similar to the lSersi"3 (|l968l ) 
law for surface densities. Differences between Eq. [TJ the 
Se rsic law and the i r stru cture parameters N can be seen 
in lTamm fc Tenied (120051 ). 

The density distributions for the visible components 
were projected along the line-of-sight, and their sum yields 
the surface brightness distribution of the model 



L{A) 



4 



li{a)a da 
(a2 - A2)V2 



(2) 



where A is the major semiaxis of the equidensity ellipse 
of the projected light distribution and Qi are their apparent 
axial ratios = cos^ Q + q^ sin^ &. The angle between the 
plane of the galaxy and the plane of the sky is denoted by Q. 
The summation index i designates four visible components. 

The model parameters q, ao, L and A' for the com- 
ponents were determined by a subsequent least-squares ap- 
proximation process. First, we made a crude estimation of 
the population parameters. The purpose of this step is to 
avoid obviously non-physical parameters - relation ((2} is 
non-linear and fitting a model to the observations is not 
a straightforward procedure. Next, a mathematically cor- 
rect solution was found for each component. Uncertainties of 
the model (related to the coupling of parameters) were esti- 
mated by using the partial s econd derivatives of the sum o f 
the least-square differences l|Bevington fc RobinsonI lioOSl ) . 
Details of the least-squares approximati on and the gen- 
eral m odelling procedure were described by|E inasto fc Ifaud| 
1 19891 ) jTenies et all (|l994l ): lTenies. Hand fc Einasto (199 j ); 
Tempel (in preparation). 

The parameters of the model (the axial ratio q, the har- 
monic mean radius ao, the structural parameters A', the 
dimensionless normalizing constants h and k, UBVRIL- 
luminosities) are given in Table [T] 



4 STELLAR COMPOSITION OF THE 
POPULATIONS 

4.1 Chemical evolution models 

The key objective of considering chemical evolution here is 
to have an estimate for M / Ls of the stellar components in- 
dependently of their kinematics. In principle, stellar M/Ls 
can be treated as free parameters in dynamical galaxy mod- 
els; masses of galactic components are then determined by 
the best fit to the kinematical data. In reality, a high degree 
of degeneracy exists in such a fitting of the components, e.g. 
the disc mass can be raised if the dark halo mass is lowered. 
Thus the use of M/Ls, provided by chemical evolution mod- 
els reduces the uncertainties arising from such degeneracy. 

In our model, galactic components are assumed to be 
homogeneous stellar populations. For this reason, in order to 
describe chemical evolution, we can use for each component 
a simple stellar population (SSP) evolution model. In SSPs, 
stars are left to chemically evolve according to the observa- 
tionally determined evolutionary tracks (colour-magnitude 
diagrams). The population is formed of born-together stars 
of various masses; the relative amount of stars of a given 
mass is determined by the given initial mass function (IMF). 

Note that unlike several other studies of both local and 
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high-redshift galaxies where a single SSP model has been 
adopted for the whole galaxy, we use a separate SSP model 
for each component, thus the galaxy as a whole is modelled 
as a superposition of several SSPs which is closer to the 
actual situation. 

A variety of SSP models have been developed elsewhere. 
We have chosen to rely on the up-to-date and much-tested 
model Starburst99, designed at the Space Telescope Sci- 
ence Institute initially for young stellar populations, but 
later adapted als o for modelling populations with ages up to 
the H ubble time (|Leitherer et al.ll 19991 : Fv^zguez fc Leithererl 
I2OO5I ). 

We constrain the metallicity of each population accord- 
ing to the measurements described in Section[2]2] We further 
constrain the models according to star f ormation history es - 
timates from the literature. Guided by lOlsen et al.l (|2006l ) 
and .Brown et al.l (2006) we suppose that the star forma- 
tion, both in the bulge and in the disc, started 10 gigayears 
ago and lasted with constant intensity for four gigayears. 
In addition we pollute the disc with a younger generation 
of stars having formed during the last four gigayears. Ac- 
cording to the observed coulours, we find the fraction of the 
younger pop u lation will be roughly 10 % . Also referring to 
iBrown et al.l ()2006t ). we set the star formation of the inner 
halo having taken place 8-12 gigayears ago and ascribe sim- 
ilar history to the outer halo. 

SSP models are also strongly dependent on the choice of 
the input IMF. While being among the most important dis- 
tribution functions in astrophysics, the IMF, on the other 
hand is rather poorly known despite the large amount of 
work devoted to it. On the basis of the field star studies 
of the Milky Way disc and star clusters, relatively narrow 
constraints on the IMF have been imposed for stars of at 
least the solar mass; much less is known about the low-mass 
end of the function. Making the matters worse, there is no 
reason to assume that all stellar populations in all galaxies 
descend from identical IMFs; birth conditions may vary sig- 
nificantly with time and environment, thus in the studies of 
other galaxies it is important to consider the effect of possi- 
ble deviations from the IMF established for our own ne igh- 
bourhood. We use the 4-segment IMF bv lKroupal l|200ll ') for 
population modelling, but also run the Monte Carlo simu- 
lations to calculate variations resulting from the provided 
IMF uncertainties. 

The effects of variations of age, metallicity and the IMF 
on the colours and mass-to-light ratios of the modelled stel- 
lar populations are demonstrated in Fig. |3] In this figure, 
the thick lines represent the evolution of a population with 
[M/H] = —0.4 (close to the average measured disc metallic- 
ity of M 31) and with Kroupa IMF; the thin lines represent a 
population with the solar metallicity (rougly the upper limit 
of the disc metallicity). The population represented with 
the dashed lines has [M/H] — —0.4, but a "bottom-heavy" 
IMF, giving rise to a higher number of low-mass stars and 
less high-mass stars (within the estimated uncertainties of 
the IMF). Fig. [3] shows that uncertainties of the IMF have 
more effect on M/L than uncertainties of metallicity deter- 
mination, and that colour indices are rather independent 
of both of them. The resultant uncertainties of the colours 
and mass-to-light ratios of the actual M 31 populations are 
demonstrated in Fig. |4]and Fig. [5] with the thickness of the 
model lines. 




2 4 6 8 10 12 14 

Time (Gyr) 



Figure 3. Evolution of B — V and I — L colour indices and 
mass-to-light ratio M/Lb at different metallicities and at differ- 
ent initial mass functions. Note that M/Lg is very sensitive to 
variations of the initial mass function. 

4.2 Model fitting 

In a number of studies, M/L or its distribution has been 
determined on the basis of just a single colour estimate. 
When the actual uncertainties of colour measurements are 
considered, the derived M/L cannot be very reliable. A sin- 
gle colour value can be fitted with a large variety of initial 
metallicity, age, IMF, etc. combinations. On the other hand, 
the use of more than one colour measurement for determin- 
ing the M/L often indicates that observations are actually 
controversial with the predicted colours and the M/L of the 
best-fitting model actually depends on which colour is used 
for fitting. Thus we have used the colour information of all 
the available observations from U to L to decrease such un- 
certainties, giving in total five independent colour-indices. 

The best-fitting model of the galaxy was found in an 
iterative process. First a model density distribution, formed 
as a sum of the four components described in Section 13.21 
was least-squares-fitted to the observed luminosity distri- 
butions. This was done simultaneously for all the available 
colours, allowing us to determine the first-order approxi- 
mation of model parameters and providing us with colour 
indices of each component. These colour indices were com- 
pared to those of the chemical evolution model for the cor- 
responding metallicity of the component, thus determining 
the first-order M/Ls for each colour. In the next iteration, 
these M/Ls were fed to the luminosity distribution model, 
fixing the ratios of luminosities in different filters for each 
component and resulting in slight corrections for the com- 
ponent parameters and also the model colour-indices, which 
were again compared to the chemical evolution models. Af- 
ter four or five iterations, the changes of the parameters 
became negligible. 

Final results of Starburst99 chemical evolution model 
fitting are presented in Fig. |4] and Fig. [5] with darker lines. 
The wide regions of a constant colour represent the observed 
colours of each population and the vertical dashed lines in 
these figures denote upper and lower limits of the acceptable 
M/L values within the age limits described in Section 14.11 
In general, the colours predicted by the chemical evolution 
model are in good concordance with the observed ones for 
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Figure 4. Integrated colours of M 31 bulge and disc components 
as a function of mass-to-light ratio M/Lg. Filled areas - con- 
straints from observations; darker areas - predictions from chem- 
ical evolution models with uncertainties; dashed vertical lines - 
limits set by population age estimates. 



Inner halo 
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Figure 5. The same as Fig. |4]but for the halo populations. Note 
that especially in the outer halo, the modelled colour stays con- 
siderably bluer than the observed one. 



all the components with the exception of the diffuse outer 
halo. Fig.[6]also shows that in general, the model colours are 
not contradicting with the observed ones, apart from some 
measurements in the outer regions along the minor axis. In 
the case of the diffuse halo, the "photometry" is actually 
based on the counts of RGB stars and is probably not very 
accurate. 



5 RESULTS AND DISCUSSION 

The absorption-corrected composite surface brightness pro- 
files in U , B, V, R, I , L colours along the major and the 
minor axes were derived by averaging the results of different 
authors and are presented in Fig. [6] with errorbars. All the 
surface brightness profiles obtained in this way belong to the 
initial data set of our model construction. With solid lines 
in Fig. [6l we show how the model fits to the restored lu- 
minosity distributions. The four-component photometrical 
best-fitting model fits all photometric profiles with a mean 
deviation - ^™°'^°') = 0.25 mag. The parameters of 

this model (the harmonic mean radius ao, the axial ratio 
q, the structural parameters N, the dimensionless normaliz- 
ing constants h and k, l7_By_R/L-luminosities) are given in 
Table [1] together with uncertainty estimates. 

Fig. [6] shows the modelled luminosity distributions 
along both the major and the minor axis. The contribution 
of the individual components to the luminosity in 1^-filter is 
presented in Fig. [7] In Fig. |8] the modelled colour distribu- 
tion is plotted. 

The total luminosity of M31 corrected for the intrinsic 
absorption is Lb = (3.3 ± 0.7) • IO^^Lq, corresponding to 
an absolute luminosity oi Mb = —20.8 ± 0.2 mag. About 
half of it (51%) is radiated by the disc, one third (34%) 
by the bulge, about 15 % by the inner halo and about one 
per cent can be ascribed to the outer halo. The dust disc 
absorbs as much as 41 % (or 0.57 mag) of the total _B-flux. 
Of the total apparent _B-luminosity, the disc, the bulge, the 
inner halo and the outer halo give 51 %, 30 %, 17 % and 2 %, 
respectively. The total intrinsic colour indices are (B — V) = 
0.9 ± 0.2, {V - R) = 0.6 ± 0.2, and {R - I) = 0.7 ± 0.3. 
Colour indices predicted partly by the chemical evolution 
model are (U — B) = 0.4 and {I — L) — 1.9. Resulting from 
the chemical evolution model mass-to-light ratios M/Lb of 
the visible components together with the corresponding ages 
and metallicities are given in Table [2] 

The total mass of the visible matter Mvis = (10- 
19)T0'^°Mq, making the mean mass-to-light ratio of the 
visible matter M/Lb = 3.1-5.8 Mq/Lq. 

The main difference from a conventional bulge -|- disc 
4- halo model is a separation of the inner, slowly rotating, 
moderately enriched inner halo and the very metal-poor, 
diffuse and extremely extended outer halo. This structure 
follows from the recent surface brightness and kinematics 
measurements of the outer parts of M31, which have re- 
vealed a slowly rotating component. The inner halo in the 
present model is actually similar to the population identified 
as th e whole halo in prev i ous stellar p opulatio n studies of 
M31 (|van den Berghlll99ll : iTenies et al.,.1994 : Geehan et al.l 
I2OO6I ). The outer diffuse halo gives only marginal contribu- 
tion to the overall luminosity of the galaxy. 

In their analysis of star counts, llbata et al.l l|2005l . l2007l ) 
derived the scale-length of the exponential "extended disc" 
along the major axis to be rd = 6.6 kpc. To match this 
value with the minor axis scale length rd = 3.22 kpc, they 
propose that the extended disc component might have an 
inclination angle of 60.8° instead of i = 77° adopted for 
the rest of the galaxy. Our model suggests the "extended 
disc" (the "inner halo" in our model) to be a thick compo- 
nent with an axial ratio of about q = 0.5, which is close to 
the value q ~ 0.55 estimated for the outer parts of M31 
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Figure 6. Luminosity distribution of M31: left-hand panel - UBVRIL (from bottom to top) luminosity along the major axis; right-hand 
panel — BVRI (from bottom to top) along the minor axis. Error bars - observations; solid lines — model. 



Table 1. Calculated parameters of the photometrical model. 



Population 


ao (kpc) 


q 


N 


k 


h 




Bulge 


0.64±0.03 


0.6±0.03 


4.2±0.3 


6.498 ■ 10-5 


4943 




Disc 


9.3±0.3 


0.05±0.03 


0.7±0.07 


0.8479 


2.240 




Inner halo 


4.0±0.4 


0.5±0.04 


2.7±0.3 


7.229 ■ 10-3 


159.8 




Outer halo 


40. 


0.9 


2.0 


0.05 


33.33 




Population 


Lu (IOIOL0) 


Lb (IOIOL0) 


Lv (IQIOLq) 


Lr (10"'Lq) 


Li (10"'Lo) 


Ll (IOIOL0) 


Bulge 


0.78±0.1 


1.09±0.1 


1.45±0.2 


1.70±0.2 


2.28±0.4 


6.62±2.4 


Disc 


1.67±0.5 


1.68±0.2 


2.02±0.4 


2.31±0.4 


3.11±0.6 


(7.88) 


Inner halo 


(0.44) 


0.50±0.2 


0.61±0.1 


0.69±0.1 


0.88it0.1 


(1.86) 


Outer halo 


(0.05) 


0.05 


0.05 


0.06 


0.06 


(0.10) 



Notes: Luminosities are corrected for intrinsic absorption and absorption in the Milky Way. Luminosities given in parentheses are not 
based on observed surface brightness profile; they are extrapolated according to the chemical evolution models. Errors for the outer halo 
are undetermined. 



by iPritchet fc van den Berghl l| 19941 ). Combined with the 
measured "warm" k inematics rather than "cold" rotation 
jGentile et aLll2007l ). we consider an "inner halo" to be a 
more descriptive name for this population. 

The total stellar mass of M 31, as estimated on t h e basi s 
of the 3.6 \im infrared luminosity by iBarmby et al. I (|2006l ). 
is 11 • 10^°Mq. This lies within our estimated range (10- 
19)-1O^°M0. A more de tailed mass d i stribu tion model of 
M 31 was constructed by[GeehanetaL| (|2006l ). who derived 
M = 3.2-lO^°M0 for the mass of the bulge, and 7.2-10"Mq 
for the mass of the disc. They also derived M/Lr = 
3.9M0/L0 for the bulge, and M/Lr = 3.3, Mq/Lq 
for the disc (not corrected for extinction). These values 
also lie within the range of our model values: bulge mass 
(4-7.5)-10^°Mq, disc mass (4.8-8.7)-lO"M0 and mass-to- 
light ratios, not corrected for extinction M/Lr = (2.4- 



4.4)Mq/Lq for the bulge and M/Lr = (1.9-3.4) Mq/Lq 
for the disc. 

Fig. |3]demonstrates that the colours predicted by chem- 
ical evolution models match well with the observed colours 
within age limits estimated in the literature. Thus we can 
conclude that chemical evolution models and age estimates 
on the basis of colour-magnitude diagrams are generally in 
very good agreement with the actual integrated colours of 
M31. However, this is not the case for the halo regions. 
Fig. [5] shows that especially for the outer halo the observed 
colours are remarkably red compared to the modelled ones; 
even after a Hubble time of evolution, the colours of the 
modelled population remain significantly bluer. What could 
be the source of this discrepancy? Quite clearly, the pres- 
ence of dust in quantities required for reddening at the level 
of E{B — V) > 0.5 mag is very unlikely in the halo re- 
gion. On the other hand, this would not be the first case 
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Table 2. Properties of stellar populations, set by observations 
and chemical evolution models. 
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Figure 7. Model surface brightness distribution in V-filter along 
the major axis and minor axis. The dashed lines indicate the 
contribution of each stellar population. 
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Figure 8. Colour distribution in M 31 along the major axis. From 
bottom to top: U — B, B — V, V — R, R — I and I ~ L. Error bars 
— observations; solid lines - model. 



Population 


Age (Gyr) 


[M/H] 


M/Lb 


Bulge 


6-10 


0.0 


3.7-6.9 


Disc 


6-10" 


-0.3 


2.9-5.2 


Inner halo 


8-12 


-0.7 


2.7-5.7 


Outer halo 


8-12 


-1.5 


1.6-2.8'' 



" Our model disc also contains 10 % of younger stars. 
M/Lb of the outer halo is based on the chemical evolution 
model only. 



of finding evidence for the presence of extremely red haloes: 
the stacked image of more than a thousand Sloan Digital 
Sky Survey disc galaxies shows sig nifican t /-colour excess 
('Zib etti. White fc Brinkmannlliooi i. while IZackrisson et al] 
(200a) have found extremely red haloes around blue com- 
pact galaxies. In the latter work, a bottom-heavy IMF was 
given as the most likely cause for the effect. 

The general ages of the bulge and disc stars suggest 
a similar formation history for these components, starting 
a few gigayears after the Big Bang and continuing over an 
extended period. The small additional younger population 
of stars, needed to achieve the somewhat bluer colours of the 
disc, are in agreement with a relatively recent merger event, 
perhaps with a structure which by now has been disrupted 
to become the Giant Stream. 

How unique is the developed model? To some extent, 
and partly because of the limitations set by the assumptions 
of homogeneous populations with axial symmetry and avoid- 
ance of ring- like structures, the components of the model are 
degenerate. With suitably chosen structural parameters, the 
spheroidal bulge and inner halo populations can compensate 
for each other. In fact, we can think of the bulge and the 
inner halo as a single component with a declining metallicity 
gradient towards outer radii. Similarly, the actual structure 
of the two halo populations is by no means fixed by our 
model, although there seems to be rather little in common 
in their nature and origin. The error estimates for our final 
model are given in Table (2] These estimates are based on 
the partial second derivatives of the sum of the least-square 
differences ( the second derivatives o f the g oodness-of-fit pa- 
rameter, see lBevington fc RobinsonI l|2003h . eq. 8.11). These 
errors describe how accurately each model parameter is de- 
termined on the basis of the used observational data. For the 
outer halo, where the observed photometry was very uncer- 
tain, these estimates were quite large (over 100 %) and we do 
not present them in Tabled The presented errors describe 
the goodness of our final model, and they do not fully take 
into account the coupling of parameters. 

We have attempted to estimate the effect of all obser- 
vational and methodological errors and uncertainties on our 
results. The dominant source for uncertainties comes from 
the estimations of the mass-to-light ratio, given by chemical 
evolution models. We have tried to combine all the available 
information on the photometry and metallicity of M31 and 
to use up-to-date evolution models of stellar populations to 
minimize these uncertainties. But at a critical glance, we are 
still left with masses and mass-to-light rations uncertain by 
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a factor of 10-30 per cent. These ranges can be tightened 
with the inclusion of kinematic information. This is done in 
Paper II. 
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